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Mean signal levels well below those predicted by standard 
techniques on low elevation earth-space microwave links have been 
observed and reported previously. Additionally, amplitude fluctu- 1 


ations under clear air conditions accompany these reduced signal 
levels and increase significantly as elevation angle decreases. 

It is shown that a simple physical model is adequate for the pre- 
diction of the long term statistics of both the reduced signal 
levels and increased peak-to-peak fluctuations. The model is based 
on conventional atmospheric turbulence theory and incorporates both 
amplitude and angle of arrival fluctuations. 


Angle of arrival statistics were compiled from previous and 
current experiments; and these results were used to establish an 
empirical relationship between the long term average angle of arrival! 
variance and either path length for terrestrial paths or elevation 
angle for earth-space paths. Secondly, the relationship between 
turbulence induced amplitude fluctuations and path length previously 
derived in the work of V. I. Tatarski and verified by measured data 
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has also been utilized. 


It is hypothesized that a plane wave propagating through the tropo- 
sphere is perturbed by the two coupled statistical processes mentioned 

above. This perturbed plane wave is then received by a finite aperture 
antenna having a Gaussian power pattern. The average realized gain is 
then calculated accounting for the angle of arrival fluctuations. The 
resulting average gain degradation, when combined with the usual 
atmospheric gas absorption, adequately predicts the reduced long term 
signal levels observed at low elevation angles on earth-space paths at 
2, 7.3 and 30 GHz. The expected value of the variance of the received 
signal amplitude fluctuations is calculated by combining the amplitude 
and angle of arrival statistics of the incident wave with the receiver 
characteristics. This model predicts the average variance of signals 
observed under clear air conditions at low elevation angles on earth- 
space paths at 2, 7.3, 20 and 30 GHz. 

Design curves based on this model for gain degradation, realizable 
gain, and amplitude fluctuation as a function of antenna aperture size, 
frequency, and either terrestrial path length or earth-space path 
elevation angle are presented. 
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The propagation of electromagnetic and acoustic waves in media 

whose constitutive parameters are random functions of time and space 
has become a topic of increasing interest. Terrestrial line of sight 
and earth-space microwave propagation paths traverse the troposphere, 
where fluctuations in temperature, pressure, and water content produce 
random variations in the electromagnetic properties. The resulting 
effects are generally more pronounced with increasing frequency and 
larger path lengths. Communications technology continually progresses 
higher in frequency and bandwidth in search of greater information rates, 
long atmospheric path lengths are necessarily employed as satellites 
come into common usage for communications relays, and remote sensing 
of meteorological, terrestrial, or airborne phenomena requires electro- 
magnetic propagation through long distances in the atmosphere. The 
design of each of these potential or current technologies requires 
knowledge of the interaction between the turbulent atmosphere and wave 
propagation based upon a realistic model of the turbulence and its 
consequences. 

The work of V. I. Tatarski [1] in 1961 appears to be the first to 
form a broad theoretical basis for plane wave propagation through a 
turbulent atmosphere. His amplitude fluctuation and power spectrum 
results were extended in 1971 to include phase and differential phase 
variation [2]. The work of Schmeltzer [3] extended this work to form 
a basis for spherical -wave, finite-aperture cases. Both sets of work 
employed the Rytov method and its inherent restrictions that the 
magnitude of fluctuations on the wave be small and that all refractive 
perturbations be large compared to wavelength. In 1969, Lee and Harp [4] 
applied a phase screen technique to this problem which permitted 
generalized amplitude and phase fluctuation expressions to include 
finite antenna apertures, focused beam waves, medium losses, and 
atmospheric anisotropy. This technique, which has a primarily physical 
rather than mathematical basis, allows the Rytov restriction of large 
perturbations relative to wavelength to be relaxed. Either or both 
methods have been variously employed by Clifford [5], Ishimaru [6], 
and Mandics, Lee, and Waterman [7] to obtain temporal frequency spectra. 
The extension of turbulence theory to account for strong fluctuations 
and explain the optical saturation phenomena was made by Gracheva and 
Gurvich [8] in 1965. 

Excellent summaries of the above work and reviews of current work 
are available [1,2,9-13]. The recent summary by A. Ishimaru [13] 
contains a particularly good bibliography of theoretical and experimental 

papers. 



In principals theoretical work in this field may be verified by 
experiment, and a great deal of effort has been made to do so. 

Experiments must be structured carefully such that all relative physical 
parameters are measured and no false assumptions are made about the 
measurement equipment or the medium. Because of the statistical nature 
of the measurement problem, limited temporal measurements of quantities 
such as amplitude, absolute, or differential phase do not adequately 
characterize random media propagation effects from the systems appli- 
cation viewpoint. Communications link design requires long term time 
characterization in the form of amplitude, phase, spectral, and coherency 
distributions. The work which follows presents a long term, time 
average, statistical model with empirical constants for microwave 
propagation through the turbulent troposphere. Amplitude and differential 
phase statistics are specifically addressed, but the techniques may be 
extended to spectral and coherency characterizations. 

The material contained in this report is also used as a dissertation 
submitted to the Department of Electrical Engineering, The Ohio State 
University as partial fulfillment for the degree Doctor of Philosophy. 
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INTRODUCTION 


The propagation of electromagnetic and acoustic waves in media 
whose constitutive parameters are random functions of time and space 

has become a topic of increasing interest. Terrestrial line of sight 
and earth-space microwave propagation paths traverse the troposphere, 
where fluctuations in temperature, pressure, and water content produce 
random variations in the electromagnetic properties. The resulting 
effects are generally more pronounced with increasing frequency and 
larger path lengths. Communications technology continually pro- 
gresses higher in frequency and bandwidth in search of greater infor- 
mation rates, long atmospheric path lengths are necessarily employed 
as satellites came into common useage for communications relays, and 
remote sensing of meteorological, terrestrial, or airborne phenomena 
requires electromagnetic propagation through long distances in the 
atmosphere. The design of each of these potential or current tech- 
nologies requires knowledge of the interaction between the turbulent 
atmosphere and wave propagation based upon a realistic model of the 
turbulence and its consequences. 

The work of V. I. Tatarski [1] in 1961 appears to be the first 
to form a broad theoretical basis for plane wave propagation through a 
turbulent atmosphere. His amplitude fluctuation and power spectrum 
results were extended in 1971 to include phase and differential phase 
variation [2], The work of Schmeltzer [3] extended this work to form 
a basis for spherical -wave, finite-aperture cases. Both sets of work 
employed the Rytov method and its inherent restrictions that the 
magnitude of fluctuations on the wave be small and that all refrac- 
tive perturbations be large compared to wavelength. In 1969, Lee 
and Harp [4] applied a phase screen technique to this problem which 
permitted generalized amplitude and phase fluctuation expressions to 
include finite antenna apertures, focused beam waves, medium losses, 
and atmospheric anisotropy. This technique, which has a primarily 
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physical rather than mathematical basis, allows the Rytov restriction 
of large perturbations relative to wavelength to be relaxed. Either 
or both methods have been variously employed by Clifford [5], 

Ishimaru [6], and Mandics, Lee, and Waterman [7] to obtain temporal 
frequency spectra. The extension of turbulence theory to account for 
strong fluctuations and explain the optical saturation phenomena was 
made by Gracheva and Gurvich [8] in 1965. 

Excellent summaries of the above work and reviews of current work 
are available [1,2,9-13]. The recent summary by A. Ishimaru [13] 
contains a particularly good bibliography of theoretical and experi- 
mental papers. 

In principal, theoretical work in this field may be verified by 
experiment, and a great deal of effort has been made to do so. 
Experiments must be structured carefully such that all relative physical 
parameters are measured and no false assumptions are made about the 
measurement equipment or the medium. Because of the statistical nature 
of the measurement problem, limited temporal measurements of quantities 
such as amplitude, absolute, or differential phase do not adequately 
characterize random media propagation effects from the systems appli- 
cation viewpoint. Communications link design requires long term time 
characterization in the form of amplitude, phase, spectral, and co- 
herency distributions. The paper which follows presents a long term, 
time average, statistical model with empirical constants for micro- 
wave propagation through the turbulent troposphere. Amplitude and 
differential phase statistics are specifically addressed, but the 
techniques may be extended to spectral and coherency characterizations. 



CHAPTER I 
METHOD 

As millimeter wave systems are applied to earth-space and 
terrestrial communications links, the designer must be aware that 

certain propagation phenomena which were negligible at lower 
frequencies will have pronounced effects on link characteristics . 

For example, a larger aperture which adds gain to the link equation 
at high elevation angle on an earth-space path may provide less gain 
than expected if employed at low angles. Similarly, clear air 
amplitude scintillations which may be neglected at decimeter wave- 
lengths on a long unobstructed terrestrial path may be quite sig- 
nificant at millimeter wavelengths. 

The temporal behavior of the signal is less important than its 
statistical properties when one is interested in link reliability. 

The long term time average behavior of signal level and variance are 
the most useful design criteria for a communication link which must 
be continuously operational over a period of years. A propagation 
model must predict this long term average behavior and provide a 
statistical measure of expected deviation from that mean in terms of 
long term probability distribution functions. 

Refractive irregularities in the troposphere due to turbulence, 
humidity, and temperature gradients affect the propagation of electro- 
magnetic waves of centimeter and shorter wavelengths quite markedly on 
long terrestrial propagation paths or on earth-space paths with small 
elevation angles. As a plane wave traverses the atmosphere, changes 
in refractive index due to turbulence, i.e., ripple or tilt, will 
perturb the phase front, induce amplitude changes by focusing effects, 
and scatter energy away from the propagation path. The following 
discussion will deal only with clear air effects and thus energy 
scattered out of the beamwidth of an antenna will be considered 
negligible compared to the energy received. Frequencies of interest 
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will range from roughly 1 to 100 GHz, which are sufficiently high that 

ionospheric effects may be neglected except in rare instances. 

The statistics of the perturbed wave are modelled by theoretical 
expressions matched to experimental observations. The perturbed wave 
impinges on a finite aperture antenna, whose response is a received 
signal with amplitude diminished from that expected for a free-space 
path and fluctuating in amplitude. 

An expression for the incident wave, the manner in which it 
attenuates and decomposes into amplitude and angle varying components 
and the transfer function of the receiving antenna will now be 
presented. 



CHAPTER II 
STATISTICAL MODEL 


An unperturbed plane wave will be considered incident on a 
turbulent medium at position z=0 (Figure 1). Propagation within 
the turbulence and the resulting perturbations of the signal 
continue until the wave impinges on a receiving aperture at Z=L. 

The incident plane wave decomposes into two uncoupled components 
as it traverses the region of turbulence. 

The first component, with magnitude f ^ , varies in amplitude 
and is constant in angle. All of the plane wave energy is contained 
in f-j at z=0, and its amplitude fluctuation is very small at the 
beginning of the turbulent region. As this component penetrates 
deeper into the turbulence, it loses energy to the second com- 
ponent and its peak-to-peak amplitude fluctuations become a larger 
percentage of its own remaining energy. 

The second component, having magnitude f^, is constant in 
amplitude at any given position z and it varies in angle or direction 
of propagation, f^ contains no energy at z=0 and at the beginning 
of the region of turbulence it is small in amplitude and has small 
fluctuations in angle. As this component penetrates into the 
turbulence it acquires energy at the expense of the amplitude 
varying component, f-j , and its peak-to-peak angle fluctuations 
increase. 

This concept of constant angle and constant amplitude com- 
ponents which are uncoupled corresponds to Ishimaru's discussion of 
coherent and incoherent waves (Reference [13], pp. 1046-1047) for 
line of sight propagation through randomly distributed particles. 

The constant angle, varying amplitude intensity <f^> corresponds to 
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Figure 1. Plane wave decomposition. 
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Ishimaru's coherent intensity, I Also, the constant amplitude, 

2 c 

varying angle intensity <f 2 > corresponds to his incoherent intensity. 


The magnitude of the electric field component, E-, of a per- 
turbed plane wave incident on a receiving aperture is a random function 
of time, varying in amplitude and phase 

^ = C[? 1 cos(w o t+C 0 +? 1 -^ 0 *r)+f 2 cos(w o t+? 0 +? 2 -e a *r)] (1) 

Random variables are denoted by tilde, and C represents the average 
magnitude of an incident plane wave, including free space path loss 
and gas absorption, after propagating through some length, L, of 
turbulent medium. f-| and f 2 are magnitude components repre- 
senting the decomposition, as a function of distance, of the original 
wave into amplitude and angle of arrival varying parts, respectively. 
The frequency of the wave is w o , is the mean absolute phase 
between the transmitter and receiver, and ^ and C, are random 
variables representing fluctuations in apparent path length. It is 
assumed that scattering out of the path does not occur and that the 
random amplitude and angle of arrival processes are conservative in 
intensity, i.e., 

<? 1 > + 4 = 1 ( 2 ) 

where < > denotes the ensemble average over amplitude and angle. 

Note that total intensity is one, that is, it is normalized to 
power density at any position along the path. This normalization 
is justified because the factor C in Equation (1) contains free 
space and gas absorption terms which equally affect f-j and f 2> In 
addition, it will later be shown that the amplitude process <?^> 
further splits into an average and fluctuating part for the case of 
strong turbulence. The mean power in the amplitude component 
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^ 2 — 2 / 2 % 2 

will take, on the form < x^ > - f-j (1+c.p , where is amplitude variance 
and f-j is the statistical mean of f ^ . This is simply a statement that, 
for the peak to peak fluctuation increases at the expense of 

power in the mean, as we saw in the sketch of f^(Z=L) in Figure 1. 

3 is the vector propagation constant, where the subscripts 0 and a 
denote the direction of propagation relative to the beam axis, r is 
the position vector. 

The receiver will be assumed linear, with ideal bandpass 
characteristics, and does not maintain absolute phase coherence. If a 
square-law device is employed as the first mixer, the receiver output 
voltage will be (see Appendix I) 

v = — ~ C |? 1 G(0)exp(-i C-, )+f 2 G(o)exp(-i? 2 ) | . (3) 

Note that absolute phase is lost and the io Q t time dependence is 
effectively integrated by the bandwidth of the receiver. G(a) is the 
directive gain of the antenna as a function of the angle of arrival, a, 
measured with respect to the beam axis. 

% % 

The first and second moments of v, namely mean <v> and second 

'V? 

moment <v >, must be obtained in order to characterize received 
signal variance and gain degradation. Hence, the statistical properties 
of f , a, L , and must be used since they comprise the expression 
for v. f-j and f^ are amplitudes of the two components and C-j and ^ 
are their respective relative phases. In the following, n and a will 
denote the random variable spaces in amplitude and angle of arrival, 
respecti vely. 


Define a as an angle of arrival random variable, referenced to 
the axis of the receiving antenna beam . Its distribution in angle 

l~2 2 

equals the mass in the circle cry\je +<j> of radius a, where e and <)> form 

an orthogonal coordinate system with the axis (see Figure 2). The 
beam axis is pointed in the mean direction of arrival of the incident 
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Of 'X* 

wavefront. Hence e and <p have zero mean angle of arrival. If they 

are also normal, independent random variables with equal variances 
2 2 2 

^ CT e = % =a 2^ 5 anc * constant amplitude dummy variable o, it may be 
shown [16] that the probability density function, pdf, of the angle 

f\j , 

varying component of v is then 

h^ ( o'- » n ) = h a (a)6(n) = ^ exp[-a 2 /2a 2 ]6(n) ; a >_ 0. (4) 

°2 

The angle varying component f,, is constant in amplitude, namely 
amplitude, and hence has zero amplitude variance. The definition of 
its pdf on the amplitude random variable space with the Dirac delta 
function <$(n) denotes this fact. The angular fluctuation in a is 
completely specified by the Rayleigh density function h (a), with 
angular variance parameter c^. Independence between the a and n 
random variable spaces is observed, since h2(a,n) is expressed as a 
simple product of pdf's, one a function of only a and the other a 
function of only n . 


The random amplitude component ?, has been found to be normally 

„ 1 p 

distributed [14]; let its mean be f ^ , its variance a-| , and express its 
randomness in terms of the amplitude dummy variable, n and the angle 
of arrival variable, a. The pdf, h^(a,n), of the amplitude varying 

component, ?-j , is then assumed to be 



expL.-- (n-f i 


2 /2a^]6( a ), 


(5) 


h j (a,n ) = h ri (n)6(a) 



10 

where 6(a) denotes that amplitude varying component f-, is consta nt 
in angle of arrival and hence independent of the angle varying com- 
ponent of v [15]. The amplitude fluctuation is expressed by the 

Gaussian density function h (n) with amplitude mean f 1 and variance 

2 11 1 % 

o-|. The assumption of a normal distribution function for f-j is based 

upon the weak scattering assumption used by Tatarskii and the fact 
that the statistical behavior of log normal and normal variables is 
identical for small fluctuations. 


Henceforth, the subscript 1 will denote statistical quantities 
in the amplitude domain n. Likewise, subscript 2 will be used to 
specify statistical parameters referred to the angle of arrival domain, 
a. This use of subscripts does not imply, however, that ?-| is only 
defined on amplitude space n, i.e., it carries constant a dependence. 

% 'X* 

The phase variables tj and ^ represent random delay variables 
caused by apparent fluctuations in path length. For reasonably long 
path lengths and assuming that all scattering mechanisms are inco- 
herent, the ^ and random delay variables will fluctuate over 
many radians and may be considered to be uniformly distributed over 
— tt/ 2 to tt/2. 


The path length dependence of f-j and f^ will now be examined. 
It has been shown [17] that the variance of the normal direction 
vector of a plane wave increases linearly with path length when 
the path length is greater than the correlation radius of the 
refractivity fluctuations. This condition is certainly met for 
paths several kilometers long. Additionally, since the axis of the 


antenna is aligned with the mean angle of arrival, the second moment 

of the zero-mean angle of arrival component follows the variance 

2 ^2 

behavior of the wavefront direction vector, i.e., a =<a >. A suitable 

2 

model for the decomposition of an unperturbed wave into incoherent angle 


of arrival fluctuating components has been suggested in Reference 13, 
Equation (2-126), to be exponential, with some characteri sti c decay rate, 
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say L , which is dependent upon the scattering cross section of the 
turbulence. The expression 

1 exp[ -L/L ] (6) 

c u 

combined with (2) will be used as a suitable model for this desired 
behavior. Also, recall that f 2 varies in angle and is constant in 
amplitude. Hence, (6) may be rewritten 

i\ = 1 - exp C-L/L q ] = f 2 (?) 

% 

All statistical quantities in the expression for v have now 
been defined and these definitions will be used to determine the 
first and second moments of the receiver output voltage, 

v = — — C | /GT^T^-| exp ( -i 5-| )+/G(&) f 2 exp(-il 2 )|, (8) 

/4tt 


If B is the beamwidth of the antenna in degrees and g(a) is the pattern 
factor, then 

and (9) 


6 (a) = g(a) (180) 2 /B 2 


* _ 1 80 A 


~= C|g(0)ft + g (cx)ff 
BvWtt, 1 d 

/ a, 

v/g(0)g(ot)f-jf 2 exp(-iC 1 +iC 2 ) 


( 10 ) 


Of 'v 


+ Vg(0)g(a)f 1 f 2 exp(+i5 1 -i£ 2 ) I 1//2 . 


Assume a representative pattern factor for a narrow beam, cir- 
cularly symmetric beam with 3 dB beamwidth, B, to be 


g(a) 


exp 


-a 4 &n 




( 11 ) 


The first moment of v is found by integrating (10) and the 
pdf h-j(a,n) and pdf h 2 (a,n) times variables containing ^ and f 2 , 

respectively. Integration on amplitude variable n is taken over 
-°° to 03 and on angle variable a from 0 to w. The phase variables 
and ? 2 are uniformly distributed over -ir/2 to tt/2 and hence 
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the integration of the last terms in (10) is nearly zero (see 
Appendix B). This integration reduces the first moment to the form 


r\j 

<v> 


1 80 A 

B/47 


c 



/g(o)n h^(n )<$ (a)dadnj 



( CO CO 

r r 


2 

+ 

v / g(a)f 9 h (a)<s(n)dadn 

L. d 



0 




1/2 


( 12 ) 


Substituting (4), (5), and (11) into (12) and using the properties 
of a Rayleigh distribution function (see Appendix B), one readily 
obtai ns 


<v> 


1 80 A 
Bf/hr" 


c 



( z \ 

2 

f 2 + f 2 
T 1 T 2 

b 


^4^n2o 2 +B 2 j 



1/2 


(13) 


Note that the mean of v approaches — — C for wide beamwidth, i.e., 

2 2 BAup 

That is, the mean received signal is not degraded since its 

beamwidth is wide compared to the turbulence effects on angle of arrival. 


The expression for v 2 is 

v 2 = 1 2° A - • C 2 |g(0)?^+g(a)f2+g(a)? 1 f2exp[-i(? 1 -5 2 )] 
B 4 7 T 


+ g(a)?-jf 2 exp[-i (l 2 -\^)\ . 


(14) 


The second moment of v is found by integrating Equation (14) and the 
pdf h^(a,n) and pdf l^o^n) times variables containing and f 2 , 
respectively. Integration on amplitude variable n is taken over 

to °° and on angle variable a from 0 to n. Again, the phase integration 
on uniformly distributed random phase variables ^ and over -tt/ 2 to 
tt/ 2 reduces the second moment to 
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a,2 

<v > 


? ? 
180 A 

B^4tt 


03 TT 

y~ /* 

[g(0)-^i h (n )6 (a)+g(a)f 2 h Q , (a)6 (n )]dadn 

0 


(15) 


It may then be shown that (Appendix C) 

( 


180 2 A 2 r 2 
<v > = A C 

B 4ir 


2 -p2 -p2 

frl+f 1+ t 2 


B 


^nZa^+B 


( 16 ) 


v 2 


A gain reduction factor R may be defined in terms of <v> and 

o 2 2 

<v> evaluated for a,-,=0 as 


R = 10 log 


<?> 2 

10 


(17a) 


a 2=0 


-2 -2 




B 


R = 10 i °g 1 0 


4«,n2a 2 +B 2 


-2 —2 
f^+f £ 
T 1 T 2 


(17b) 


This gain reduction represents available received mean signal power 
relative to the signal power available if angle of arrival effects 
were absent. A similar measure of received signal variance expressed 
in dB is 


= 10 log 


10 


0,2 o 2 
<v >-<v> 

0 2 

<V> 


10 log 1Q 


—? 2 
r k 

-2 2 2 

a- x 

1 1 8^n2a 2 +E? 


! ^4%i2a 2 +B 2 


-2 —2 

T + f ^ 
T 1 r 2 


B‘ 


4«-n2a(;+B 


(18) 


(19) 



This signal variance is a measure of fluctuating power in the received 
signal compared to the available received DC signal power. 


The two statistical quantities which have now been defined, namely 

2 

gain reduction, R, and signal variance, S , may be incorporated into 
distribution functions which are of the form used in link design. The 
long term time behavior of the received signal on a microwave link may 
be characterized by a fade distribution function. The turbulence 
induced effects with which we are dealing are long term average 
phenomena and hence are present at all times. A hypothetical low 
elevation angle fade distribution is presented in Figure 3. The abscissa 
is referenced to the signal level received in the absence of turbulence, 
i.e., including free space loss and gaseous absorption. The point at 
which the signal level is R dB is also the mean of the received signal; 
thus, one point on the fade distribution is established. The fade 
distribution for turbulence induced fluctuation is assumed to be log- 
normal, with mean and median being equal. The fade distributions 
resulting from the Ohio State University ATS-6 30 GHz beacon measure- 
ments (Ref. [30], pp. 72-75) indicate that this log-normal assumption is 
valid for elevation angles above approximately 2°. A similar observation 
was made concerning the 7.3 GHz fade distributions above 4° elevation 
angle observed by McCormick and Maynard (Ref. 35). 


A fade distribution may now be produced using this assumption of 
linearity. Referring to Figure 3, it was noted that the point at which 
the received signal level is R dB represents the mean signal level. For 
a normal distribution, the mean is plotted at the 50% time abscissa 
exceeded point, indicated by (T) in Figure 3. One standard deviation 
to the right of the mean on a normal distribution occurs at the 15.9% 
time abscissa exceeded level. It may be shown (Appendix E) that one 
standard deviation of received signal level, expressed in dB and 

O 

denoted o v ^, may be written in terms of the signal variance S . This 
point, 0 v c jb bo the right of R is denoted by ( 2 ) in Figure 3. A 
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Figure 3. Fade distribution function. 
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straight line drawn between points CD and (?) now approximately 
represents the fade distribution, referenced to the mean signal level 
in the absence of turbulence induced fluctuation. This distribution 
was based on small fluctuation arguments and should be employed as a 
lower bound when estimating a particular fade distribution. 

Deviation of this fade distribution from the expected linear 
form will occur at small time percentages. Additional fading due to 
precipitation, abnormal refraction, or inversions in the atmosphere 
will cause greater fade depths for the small time percentages. However, 
the turbulence effects, which are always present, are still dominant 
for larger time percentages. For high elevation angles, i.e., short 
path lengths, S will be very small and the line drawn through points 
CD and (2) will be virtually vertical. However, the precipitation 
effects at the lower time percentages will still be present for short 
path length cases and will become the dominant feature of the fade 
distribution. 

O 

The expressions for R and S will eventually be compared with 
measured data as a function of path length. But first, it is necessary 

p p 

to obtain expressions for and 



CHAPTER III 
EMPIRICAL CONSTANTS 


Several measurements of angle of arrival variability of electro- 

magnetic waves as they propagate through a turbulent atmosphere have 
been made. Ten experimental papers dealing with these measurements 
were examined and compared [18-27]. All results are for terrestrial 
microwave links operating at frequencies from 1 to 35 GHz over paths 
ranging in length from 5.5 to 80 km and located in various climatic 
regimes. In addition, angle of arrival statistics for an earth-space 
path from the Ohio State University CTS (Communications Technology 
Satellite) 11.7 GHz Beacon Experiment [28,29] are included. 

Comparison between terrestrial and earth-space path propagation 
data may be made only if one assumes an equivalent atmospheric path 
length for the earth-space link. Although vertical gradients of 
temperature, pressure, water vapor density, and hence, refractivity 
exist in the atmosphere, one may replace this real atmosphere with a 
homogeneous atmosphere of some given height under standard temperature 
and pressure conditions for long term statistical purposes. It has 
been found, using amplitude scintillation measurements of the ATS-6 
(Applications Technology Satellite - 6) 30 GHz beacon, than an equiva- 
lent atmospheric height of 6 km is sufficient for the prediction of long 
term time average scintillation behavior. Under this assumption, the 
CTS data taken at a 32° elevation angle traverses an atmospheric slant 
path over a curved earth equivalent to a 11.3 km terrestrial path. 

The statistic to be compared from the several sources of data is 

the variance of the angle of arrival with respect to the mean direction 
2 

of arrival, a ^ expressed in dB below square degrees as 

a? = 10 loq in o? . (20) 

dB 

The normalization with respect to one square degree is arbitrary but is 
necessitated by the fact that normalization to the mean angle of arrival 

has no physical significance. The papers in which these data appear 
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are not consistent in their definitions of variability. Hence, the 
following three conventions were assumed in an effort to make a meaning- 
ful comparison of the results 


1) 


2 ) 


if variability was expressed as standard deviation a ^ in 
parts per million (ppm) of wavelength (a), 


- 10 


-dB 


iog 10 


Arctangent 


Og x A 


( 21 ) 


iver Separation 
if variability was expressed as a cumulative distribution 
function of relative angle of arrival, a , and plotted on 
Gaussian probability paper, 

at 2.5 percentile 


2 

0 2 


dB 


= 10 log 


10 


a at 97.5 percentile 

( 22 ) 


3) if variability was expressed as the observed maximum peak-to- 
peak excursion in a. 


Oo 10 log in 
aB IU 


0.95 x max. peak-to-peak excursion of a 
4 

(23) 


2 

Table I summarizes the data and the method used to arrive at c^. The 

resulting angle of arrival variances are also shown in Figure 4. 0b- 

o 

taining o ^ is quite subjective for many of the data sources, especially 
where limited measurements were presented. When a range of variances or 
probability density standard deviation limits were available, those limits 
were plotted as vertical lines through the data points in Figure 4. 


The data of Table I suggest that variability increases proportion- 
ately with path length. Hence, the angle of arrival variance in dB 
ag was plotted versus path length L in Figure 4 on a logarithmic scale. 
A first order regressive fit was then performed, resulting in: 

a 2 = 7.07 x 10" 6 x L 1 ' 56 (deg 2 ) (24) 

with L expressed in kilometers. This relationship is shown as a 

sol id line in Figure 3. 
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TABLE I 


Fig. 4 
Symbol 

Data 

Source 

Freq. 

(GHz) 

Path Length 
(km) 

2 

op Average 
(aB below deg^) 

Method 
(1 ,2,3) 

• 

Deam & Fannin 

9.35 

5.5 

-38.2 

1 

o 

Herbstreit & Thompson 

1.05 

6.5 

-37.0 

3 

□ 

Lees 

35 

10.5 

-36.5 

3 

A 

O.S.U. 

(equivalent 
Earth-Space path) 

11.7 

11.3 

-32.8 

2 

t 

Lai-iun Lo 

15.3 

15.2 

-40.75 

3 

m 

Etcheverry et.al . 

93 

18.95 

-31.4 

3 

0 

Lee & Waterman 

35 

28 

-25.2 

3 

+ 

Akiyama et.al. 

11,24 

29 

-27.7 

2 

+ 

Akiyama et.al. 

11,24 

33 

-27.4 

2 


Bell 

11.0 

55.4 

-32 

2 

X 

Janes et.al. 

9.6 

64.25 

-20.4 

1 

X 

Janes et.al. 

34.5 

64.25 

-22.9 

1 

+ 

Akiyama et.al. 

24 

78.7 

-21.62 

2 


Funakawa et.al. 

12.6 

80 

-21 .6 

2 



on 



PATH LENGTH ( KM ) 


Figure 4. Measured angle of arrival data. 
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Theoretically, it can be easily shown [17] that angle of arrival 

variance should increase linearly with increasing path length in a 
turbulent atmosphere under the assumption that L>>l_ o , i.e,, the path 

length is much greater than the scale size of the largest turbulence- 

induced inhomogeneities present on the path. Also, angle of arrival 
2 

variance is mainly due to the large-scale components of the 
turbulence spectrum. Under these assumptions the turbulence theory 
of Tatarski leads to the following expression for 

a\ = .0885 C 2 d“ 1/3 L (deg 2 ) (25) 

where C is the refractive index structure constant. This theoretical 
n 

expression is also plotted on Figure 4 for representative values of 

C 2 , where d is aperture diameter and was taken as 1 m for the theoretical 

curve. The form of Equation (25) assumes that d<</AL, a valid assumption 

for millimeter and longer wavelengths at path lengths of several 

2 

kilometers. Several uncertainties in the parameters of a 2 may account 
for differences between the theoretical and the empirical regression 

o 

curves. The angle of arrival variance a 2 may vary ±1.3 dB due to 
differences in aperture diameter d, ±18.2 dB due to typical variations 
in the atmospheric structure constant C , or ±4.6 dB due to the range 
of path elevations encountered in the experiments [31]. Additional 
differences may be due to inappropriateness of the plane wave assumption, 
especially for short path lengthsorthe requirement that L>>L o , again, 
questionable for short path lengths. 

The earth-space slant path data of Ohio State University falls 
within 2 dB of the regression line and, within the same limitations as 
the terrestrial data, agrees well with both the theoretical result as 
well as the other terrestrial measurements. 

The empirical relationship (Equation (24)) which has just been 

presented obviates the problem of calculating a 2 for our model in terms 
of a "representative" value for C 2 , which may range over two orders of 

magnitude from time to time. 
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2 

The last parameter to be determined is a - |, the variance of tSie 

amplitude fluctuations. From Tatarskii [32], one finds that for weak 
fl uctuations 


a 


2 

1 



k 7/12 l 11/6 


(26) 


where K is the propagation constant. 

2 2 

We now have expressions for a- j and a^, both having unknown con- 
stant multipliers, and an unknown length, L , each of which must be 

empirically determined in order to predict long term average behavior 
2 

of variance S and gain reduction R. In the next section measured 
variance data were used to determine these constants. 

2 

Variance S was measured as a function of elevation angle at 

2 GHz with a 10 meter diameter antenna and at 30 GHz with a 5 meter 

diameter antenna [30]. A regressive power law fit was available for 

the data from elevation angles of 2° to 40° for the 2 GHz data and 

.5° to 40° for the 30 GHz data. These limits are such that the horizon 

did not lie within the 3 dB beamwidths of the antennas. Expressions 
2 2 

for 0 -j and in terms of the three unknown constants were used in the 
model of (Equation (19)). This expression was compared to the 
regression line from measured data at 2 and 30 GHz and the three con- 
stants were varied to obtain a regressive fit minimizing the following 
mean-square-error: 
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measured 
30 GHz 


modet* 
30 GHz 


1 

el evati on 

_ 4 v nco 

- i X . uu 



200 



1 

E t 2 - 
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i=5 measured 

model 



2 GHz 

2 GHz 

✓ 

_ "J ~ , 


= ix . 05 c 


(27) 


Expressions (24) and (25) with both |J*^ an d |J'^ dependence for a 2 

were employed in the regression defined above. The minimum error 

obtained for the L ‘ dependence was approximately 40% less than that 

for the lJ ,( ^ dependence case. Hence, although theory predicts linear 

behavior of angle of arrival variance as a function of patli length, the 
1 56 

L ‘ form was used to compare the model to measured data and to obtain 
design curves. The final constants derived by the procedure outlined 
above which are used in the remainder of this work are: 


L q = 180 (km) 


(28a) 

o 2 = 2.6 x 10- 7 F 7712 


(28b) 

o 2 = 5.67 x IQ' 6 L 1 ' 56 <r 173 

(deg 2 ) . 

(28c) 


Now that the model is complete, results from it will be compared with 
measured gain degradation and amplitude scintillation results in the 
next two sections. 



CHAPTER IV 
AMP LI i UDE VARIANCE 

The expression for received signal variance (Equation (19)) was 

compared as a function of elevation angle with data acquired from the 
20 and 30 GHz ATS-6 satellite experiment as the satellite was moved 
from equatorial orbit at 94°W longitude toward 35°E longitude in 1975. 
The ground terminal at The Ohio State University, Columbus, Ohio 
employed a 4.5 m diameter parabolic, shared aperture antenna for both 
20 and 30 GHz. The linear receiver output voltage was sampled at a 
data rate of 10 samples/sec and the variance, expressed in dB with 
respect to the DC power level for N samples, was calculated from 


= 10 log 


10 


I (v. - v)‘ 
i = l 1 

-2 
N \T 


where 


N 

I v. 
1 = 1 1 
N 


n - 1024 (29) 


(30) 


The apparent elevation angle of the satellite varied from 42° to 
almost 0° as it drifted eastward, and there were no optical obstructions 
along the propagation path. Time average variances, calculated according 
to Equation (29), were recorded for each elevation angle, consisting of 
clear air data periods of one to several hours in duration. Atmospheric 
path lengths were calculated assuming a spherical earth with a homo- 
geneous atmosphere of height, h, and effective radius, R, taken as 
4/3 the actual radius, or 8479 km, to account for standard refraction. 

The path length is then 

L = [h 2 + 2hR + R 2 sin 2 e] 1/2 - R sin e (31) 

where 0 is the elevation angle. 
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The data were fit in a minimum mean-square error sense to the form 


S 2 - AL B (32) 


using expression (31) for L and also regressively fit to determine 
the effective height, h, The results were 


s 2 

b 20 

= 10 " 6 - 39 L 2.29±.l 

(33) 

S 2 

j0 

= 10 - 6 - 2 L 2.35±.l 

(34) 

and h = 

5.9±1 km . 

(35) 


The model for received signal variance (Equation (19)) was then compared 
to these results using path lengths also calculated from Equation (31). 
This comparison of the calculated results and the regression fit to the 
experimental data is shown in Figure 5. Notice that the agreement is 
within a decibel or two over most of the range shown and for very low 
elevations the magnitude of the 30 GHz fluctuations is approaching the 
level of the DC component of the received power. 

A second set of data was acquired in 197G as the ATS-6 satellite 
was returned to 94°W. A 2 GHz linear receiver utilizing a 9.1 m 
diameter parabolic antenna was implemented and variance data were 
obtained according to Equations (29) and (30). The 20 GHz transmitter 
was not functioning during the 1976 transition, but the 30 GHz beacon 
was available. The comparison between the amplitude variance 
theoretical model and the 2 and 30 GHz experimental data is shown in 
Figure 6. The agreement between the theoretical model and the experi- 
mental results is better than that shown for the 1975 data. 
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A final example compares the theoretical variance model with 
measurements of the IDCSP x-band beacon reported by R. K. Crane [33]. 

This experiment employed the 60 foot parabolic antenna at Westford, 

Massachusetts operating at a frequency of 7.3 GHz. Received power 
fluctuations were observed on angles from 0.5 to 10.0 degrees. Maximum 
and minimum limits were estimated from the data points presented in 
Reference [33], Figures 19 and 20, and plotted consistent with the 
definition of S . The mean of the limits, in dB, was also calculated. 
Figure 7 presents the means and limits of the measured data along with 
the variance predicted by the theoretical model. Again, the agreement 
is quite good. 

The model (Equation (19)) for variance consists of amplitude and 
angle of arrival induced components, with each mechanism having greater 
or lesser importance depending on the path length, frequency, and 
beamwidth. To illustrate, consider Figure 8. The model for a 30 GHz 
signal and 4.6 in parabolic antenna is presented along with the constituent 

o 

variance components due to amplitude only (i.e., c 9 = 0 in Equation (19)) 
and angle of arrival only (i.e., = 0 in Equation (19)). Note that 

the amplitude mechanism is most important at higher elevation angles and 
increases as elevation approaches zero. However, at about 11° the 
angle of arrival mechanism equals that of amplitude and then becomes 
dominant as the elevation angle decreases toward zero. Hence, amplitude 
scintillation resulting from angle of arrival effects for this case are 
relatively insignificant above 11°. 

The manner in which the incident wave splits into angle and amplitude 
components is illustrated in Figure 9. The case plotted is for a 
4.6 m diameter antenna at 30 GHz. The intensity has been normalized 
to the total power received at any position along the propagation path, 
i.e., atmospheric gas loss and free space path loss are not included. 
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Figure 7. IDCSP X-band measured amplitude variance 
compared to theoretical model. 



Figure 8. Amplitude and angle of arrival components 
of amplitude variance. 


F * 30 GHz 
d s 4.6m 



PATH LENGTH (KM) 

Figure 9. Amplitude and angle of arrival constituents 
of incident intensity. 

2 

The amplitude intensity, <f,>, dominates for path lengths less 

2 

than about 120 km in this case. For very long paths, <f-j> becomes rather 

small and, in addition, power begins to transfer into the variance 
2 2 2 

component a-|<f-j>. The angle component <f2 > dominates for path lengths 
longer than 120 km thus implying the dominance of angle of arrival 
effects on received signal variance at low elevation angles. 

The comparisons between the ATS-6 and IDCSP measurements and 

the theoretically predicted amplitude variances tend to justify the 

2 

empirical results used to model amplitude variance a, and angle of 

2 

arrival variance a^. Although the method used to obtain average angle 
of arrival statistics as a function of path length involved subjectively 
reducing numerous reported forms of data to a single definition of 



32 


2 

o' the resulting theoretical model seems to predict long term time 

L. 

average amplitude fluctuation quite well. Daily variations of the 

atmospheric structure constant are averaged out in this type of 

modeling. However, the system designer may wish to account for the 

2 

response of his communication link to the full range of expected C 
variation which has been observed to be nominally ± 10 dB about the 
long term mean [34], An example of such a maximum, minimum, and average 
received signal variance is presented in Figure 10, again, for a 4.6 m 
diameter antenna operating at 30 GHz. The maximum ± 10 dB range is 
present at higher elevation angles but shifts somewhat lower with respect 
to the average as the elevation angle decreases to zero. Similar behavior 
is expected for other combinations of frequency and aperture size. 

The frequency dependence of the amplitude variance for a fixed 
aperture size is shown in Figure 11. This family of design curves 
presents long term time average received signal variances for a 4.6 m 
diameter parabola with frequency ranging from 1 to 100 GHz. The com- 
parisons between theoretical and experimental results presented above 
indicate that the model is acceptable for frequencies betv/een 2 and 30 
GHz. The curve for 100 GHz is included only as an indication of the 
degree of amplitude scintillation to be expected if extrapolation 
using this model is warranted. This model establishes a lower bound 
on the expected amplitude scintillation at 100 GHz because refractive 
layers, abnormal focusing, or non-negl igible scattering at this fre- 
quency, and consequently very narrow beamwidth (.04°) for a 4.6 m 
aperture, will enhance the scintillation. 

Comparison of measured and calculated variance lends a measure of 
credibility to the model, at least from 2 to 30 GHz, The second quantity 
which the theoretical model predicts, namely gain degradation, will now 
be compared with experimental data as a function of elevation angle. 
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Figure 10. Maximum and minimum effects of C on 
amplitude variance. n 
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Figure 11. Frequency dependence of 
amplitude variance. 


CHAPTER V 
GAIN DEGRADATION 


The expression for long term time average received signal level 
(Equation (17)) as a function of elevation angle was also compared with 
measured satellite data. As the ATS-6 satellite (2 and 30 GHz beacons) 
was moved from equatorial orbit at 35°E longitude toward 94°W longitude, 
average received signal level was recorded and plotted as a function 
of elevation angle from 0° to 40° [30]. A 4.6 m diameter Cassegraini an 
antenna was employed at 30 GHz and a 9.1 m diameter focal point feed 
antenna was utilized at 2 GHz. In addition, median signal level as a 
function of elevation angle is available from measurements made by 
McCormick and Maynard at the Communications Research Center in Ottawa, 
Canada using the US TACSATC0M-1 7.3 GHz beacon [35]. The median signal 
level was received with a .3° beamwidth antenna as the satellite 
drifted westward, with elevation angle decreasing from 6° to 0.5° over 
a period of 23 days. The data from Reference [35], Figures 1 and 2, 
were presented as a series of distributions of received signal level 
as a function of one degree increments of elevation angle. The means 
were assumed to be the signal levels at the 50% time ordinate and 
these means were associated with the elevation angles at the center 
of the one degree increments. 

The theoretical model was used to predict mean signal level 
degradation due to atmospheric fluctuation as a function of elevation. 
This gain degradation was combined with atmospheric gas loss calculated 
for a 6 km equivalent height homogeneous atmosphere at standard 
temperatuare and pressure and plotted in dB relative to the received 
signal level at 90° elevation angle. Predicted curves are presented 
in Figure 12 for 2, 7.3 and 30 GHz for antenna beamwidths of 1.8°, 0.3° 
and 0.15°, respectively, along with measured mean signal levels from the 
ATS-6 and TACSATC0M experiments. The agreement is quite good, verifying 
that the theoretical model adequately predicts long term time average 
gain degradation as well as amplitude variance. 
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An additional plot, Figure 13, presents only the gain degradation 
component due to atmospheric turbul ence , i.e., excluding atmospheric 
gas loss. Compared to Figure 12, gas loss is certainly the major contri- 
bution to signal loss at low elevation angles, but gain degradation 
due to atmospheric turbulence is also significant, especially for the 
narrow beamwidth , millimeter wave-length case. It may be shown that 
the most significant parameter in the fluctuation degradation component 
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Figure 13. Gain degradation component of 
received signal level. 

is beamwidth. To illustrate, consider Figure 14. Gain degradation, due 
to atmospheric fluctuation only, as a function of elevation angle is 
presented for several beamwidths at a frequency of 30 GHz. Less than 
1 dB of degradation occurs down to elevation angles of about 10° 

(34.2 km equivalent path length) even for 0.05° beamwidth. Hence, the 
long term time average gain degradation is relatively minor for path 
lengths less than 34 km or elevation angles above 10°. To facilitate 
an estimation of expected gain degradation in communication link 
design, several design curves will now be presented. 
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Figure 14. Beamwidth dependence of gain degradation. 


Realized gain as a function of antenna beamwidth or equivalent 
aperture diameter at 30 GHz is plotted in Figure 15. All equivalent 
aperture diameters, given a particular beamwidth and frequency, will 
be presented for an antenna efficiency of 0.6. The curve representing 
zero path length L is simply the common gain approximation 


G = 


41253 


( 36 ) 


Realized gain curves for path lengths of 50 to 300 km are plotted 
using the theoretical model. Equivalent earth-space path elevation 

angles assuming a 6.0 km height homogeneous atmosphere are presented 
in parenthesis. 
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Figure 15. Path length dependence of gain degradation. 
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Notice that gain degradation due to turbulence induced fluctu- 
ation is negligible for beamwidths wider than about 0.7° for all path 
lengths. Degradation effects then gradually increase as beamwidth 

narrows from 0.7° to 0.05° and at any particular beamwidth are 

approximately directly proportional, in dB, to path length. As beam- 

width narrows beyond 0.05°, a saturation effect occurs and the 

degradation becomes constant for any one path length. An elementary 

explanation is that for any path length, as beamwidth becomes extremely 

small, the angle of arrival (incoherent) component of the incident 
2 

wavefront, <f£>, contributes little to the power received by the aperture 
since this angle of arrival component has negligible probability of 
being within the 3 dB beamwidth of the antenna. However, the amplitude 
component (coherent and on-axis) remains constant and the power which 
the aperture receives due to this component is only a function of path 
length and not aperture size. 

The frequency dependence of gain degradation is illustrated for 
several beamwidths at 3 and 30 GHz in Figure 16. For the longest path 
length, 300 km and for beamwidths wider than 0.2°, gain degradation is 
virtually independent of frequency. At the narrowest beamwidth, .05°, 
the degradation differs by approximately 1 dB. Hence, the turbulence 
induced degradation effect is quite insensitive to frequency, as 
opposed to the strong dependence it has on beamwidth. 

The data of Figure 16 may be quite useful to the design engineer 
employing large apertures at millimeter wavelengths. For example, 
after considering atmospheric gas loss at a particular frequency on a 
low elevation earth-space path for a 0.2° beamwidth antenna, the 
engineer may decide that he must double the aperture diameter (halve 
beamwidth) in order to attain a particular system margin. However, 
looking at Figure 16, if the equivalent path length were 200 km the 
link design should include a 4 dB average gain loss at .2° beamwidth 
and still another 4 dB when going to .1°. In some cases, it might be 
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Figure 16. Frequency dependence of gain degradation. 

more economical to employ self-phased array antennas having elements of 
relatively wide beamwidth (say .5°, referring to Figure 16) such that 
adequate gain is achieved by adding elements whose individual power 
patterns are wide enough so that average gain degradation and amplitude 
variations due to turbulence are negligible. 


CHAPTER VI 
CONCLUSION 


A model including the two coupled mechanisms of microwave angle 

of arrival fluctuation and amplitude scintillation due to tropospheric 
turbulence appears to be adequate for the prediction of long term time 
average received signal levels and amplitude fluctuations on earth- 
space and terrestrial propagation paths. The utility of an equivalent 
6 km high homogeneous atmosphere and a long term expected value for 
the atmospheric structure constant has been emphasized when deriving 
the empirical constants necessary to produce design curves from the 
model. The model does not address orographic or marine effects, but 
its ability to predict observed data from Massachusetts, Ohio and Ontario 
indicates that long term statistics of atmospheric turbulence may 
not be strongly dependent upon the particular location or climatic 
regime. 

Design curves based on this model indicate that both received 
signal level reduction and amplitude scintillation due to tropospheric 
turbulence are most apparent at narrow beamwidths. In the design of a 
long path length terrestrial microwave link or a low elevation angle 
earth-space link, the engineer must consider gain loss and amplitude 
scintillation due to tropospheric turbulence if he wishes to use large 
aperture, narrow beamwidth antennas. In cases where the angle of 
arrival mechanism dominates, it may be advantageous to utilize self- 
phased arrays to circumvent the gain degradation and scintillation 
introduced by this mechanism. The theoretical model, when properly 
applied, will indicate when such measures are necessary. 
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APPENDIX A 
LINK PARAMETERS 


The magnitude of the electric field component |E^ | incident on an 
antenna aperture is given in the text in Equation (1). 

The power density at the receiver, P , in terms of the medium 
impedance, Z , may be expressed as 

p = 1 ^ n 
r 2 L ' o 

where C is the magnitude of the incident electric field intensity. 

The following parameters are required to characterize the com- 
munications link: 


1. transmitted power in watts 

2. G t transmit antenna power gain 

2 2 2 

3. free space loss = 16rr L /A for a path length L 
at wavelength A 

4. L atmospheric gas loss (water vapor & oxygen) 

2 

A g = A G(a)/4ir is the effective aperture area of the receiving 
antenna with parameters 

1. A wavelength 

2. G (a) antenna power gain (a function of angle a). 
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Alternatively,, the antenna gain may be expressed in terms of a pattern 
factor, g(a), defined by: 

G (a) = — X -5-^2— g (a) , 
tt Ep 

where B is the antenna half-power beamwidth expressed in degrees. A 
representative Gaussian antenna pattern factor (Equation (11)) may be 
expressed as 

g(a) = exp 

g dB ( “) = - 3 

Gain degradation, R, as discussed in the text, results from angle of 
arrival fluctuation and is a function of frequency, path length, and 
aperture size. 


cc 2 4&n2 


(B) 

L_ ' ? ' _J 


All of the parameters defined above may be written concisely in 
terms of decibels to describe received power in terms of transmitted 
power, transmit antenna gain, free space loss, atmospheric gas loss, 
gain degradation, and effective aperture area as: 


P 


r 


dBW 


t 

dBW 



U - L . 
fs atm 

dB dB 


R dB + A e 


dB 



APPENDIX B 
FIRST MOMENT OF v 


The receiver output voltage may be written: 


v = — C |/g(oT ^expl-i^) + /G(a) f 2 exp(-iL) | , 

v 4 tt 


V = 180A 




C|/g(oT?iexp(-i^i) + / g( a) f 2 exp( -i e 2 ) | . 


*\j % *\j 


v is thus a function of random variables a, n, ^ , and ? 2 , or 


y 180 A r ^ y % 

V = C 0( a, n 


Bv^Ttt" 


CNaV.^yi . 


= f^To) ^expH^) + /gW f 2 exp(-i£ 2 ) 


The expected value of v is found by taking the ensemble average over 
random variable spaces a, n, ^ , and ? 2 : 

v = —z C<|n(<x,n,5,,S 9 |> 

B/4V 1 ^ 


'b 

The form of this expression for <v> is not integrable in general. 
However, two useful properties of the kernel, )o| 5 will enable an 
approximation to be invoked which results in an integrable form. First, 
the magnitude of the kernel |o| is large compared to the size of the 
expected fluctuations of ft. This condition is due to the fact that the 
development of the model was based upon the weak scattering assumption. 
Second, the kernel |n| is slowly varying over the domains of a, n, , 
and 52* For example, with n, 5], and fixed, let a vary over its 
range. The fluctuation which a. imposes on ft as it varies, and hence 
the fluctuation of |fl|, will be proportional to ] /g ( ocY | , a smooth 
function across the fluctuation domain. Similar observations may be 
made as |ft varies in o as |f ] | , in as |exp(-i 5] ) I , and in Ki as 

|exp(-i? 2 ) . 
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It may be shown (Reference [36]) that the first moment of a 
function of a random variable may be approximated by the function of 
the first moment of that random variable under the two conditions 
mentioned above. In this case & is a random variable which is a 
function itself of independent random variables n, a, £] , and 52 and 
has a well defined moments. Assuming that the small fluctuation and 
smoothness conditions are met and noting that the function of 0 is 
absolute value, then 

< |n|> = <^[m*]^ 2 ^> 


Assuming this to be the case in the expression for 


% 

<v> 



180A 

<V> 

B/4t rr 


1 80A 

<v> 

B/4tT 


1 80A 

<v> 

B/4tt” 


C<|«(a,n,C 15 5 2 ) |> 


C<|/g(o) ^exp(-i ^ ) + /g(< “T f 2 exp(-i ? 2 ) ! , 

C<[ (v / gTo)f -j ) 2 + (/g‘(a)f 2 ) 2 + ? 1 f 2 exp(-i? 1 +U 2 )g(a) 


'X, 


+g(a) f 1 f r 2 exp(+i ? 1 -i £ 2 ) 3 1 /, 2 J> 


the ensemble average may be taken as very nearly equal to the averages 
of the constituents of the absolute value process - 


<v> 


1 80 A 
b/4tt" 


Cl^gToyf ^ 2 + < C» / g(a)f z y Z + <^ 1 f 2 exp(-ic 1 +i? 2 )g(a)> 


<g(a)^ 1 f 2 e xp(+i? 1 -it 2 )>] 1/2 . 

The ensemble averages include integration on 5 ] and 52 taken over the 
domains -tt/2 to ir/2. Since and are uniformly distributed over 
this range, the last two ensembles are zero. Integration on the 
remaining two terms involves a over the domain 0 to 7 r, and integration 
on n from - M to «. Explicitly, 
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<\j 

<v> 


180 A 
b/4tt 


\2 


' / gl0)n ( r, ) 6 (a) dadn 


V -« o 



TT \2 

+ 1 I I /g(a)f^6 (n)h (a)dadri 


a 


o 


1/2 


The integration of the delta functions gives: 


<v> 


1 80A 


2 /i 


11/2 


^g(o) n h n (n)dnj ■+] 


/g(a)f2h a (a)da 


Using the fact that g(a=0)=l and f^ is independent of a and has mean 
f^, substitute the expressions for g(a), h , and h^: 


0, 

<V> 


180 A 

B/4tc 



r ^ 2 

\ n exp[-(n-f ) 2 /2o 2 ]dn 

/ 


, 1/2 


f 

TT 

2 1 


2 

f 

f 2 ’ 

exp 

a 4*n2 
2 

a 

—y exp - 

a 

„ 2 

da 

v - 

0 

2l> 

°2 

U 0 zJ 

L 


Random variable has been replaced by its dummy amplitude random 
variable, that is n. The first integral is the mean of a Gaussian 
random variable ^ , i.e. f ^ , and the second integral may be rewritten 
so that: 


a, 

<v> 


1 80A 

b/47 


f? + 


f 



(4&n2 2 +B 2 )a 2 

1 

2-j 


L. 

_ 2„2 

da 


°2 

2a ^B 

1 



1/2 
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We have assumed that «<<1 radian, so there will be negligible contri- 
bution to the integral on a over the range 7T<a<“. Therefore, let the 
upper limit of integration go to infinity and define: 

1 Un2o ? 2 +B 2 

~~~o ~ 

2r 


For a Rayleigh distribution, 


oo 



a exp[-a^/2A^]da 


A 2 - \ a exp[-a^/2A^]da. 


Rewrite the integral in the expression for <v> in terms of A: 


< v \ = M. c 

B/ 4ir 


T i +I 


h 


-I- \ a exp[-a^/2A^]da 


V 


\ a o -L 
V 2 o 


) 


1/2 


One immediately notes that the integral may be replaced by A : 

ll/2 

I / tW 

<}> = IMA, c 



L ^ 
f 2 2 

2 


°2 ) 



Substituting in the definition for A , the first moment of v is 


<v"> = c 

bAV 


/ B 

'l + '2 2 2 

1 ^W^c^+B 


1/2 




APPENDIX C 

SECOND MOMENT OF v 


With the receiver output voltage equal to 


v = 180x 


B/4tt 


C^Ta) 


[? 1 exp(-i? ] )+f 2 exp(-i ? 2 )]- 
[ ^1 exp (+i g-| )+f 2 exp(+i? 2 )] 


1/2 


v may be written as: 

v 2 = C 2 g(a) |?^+f 2 +?- ) f 2 exp(-i C 1 +i? 2 )+? 1 f 2 exp(+ii' 1 -i? 2 ) | . 

B 4ir 


*\> f \jO 

The second moment of v, namely the ensemble average of v taken over 
random variables n, a, and 5 is: 

< v 2 > = ] 2 ° A C 2 <(g(a)[? 2 +f 2 +? 1 f 2 exp(-i? 1 +i5 2 )+^f 2 exp(+i5 1 -i? 
B 4tt 

. y 'Xj 'x, 

Since fp f 2 , Cp and ? 2 are statistically independent, we may write: 

< v 2 > = C 2 [<J(0)? 2 V<J(a)f2V<Ji(a)?' 1 f 2 exp(-ic 1 +ic 2 )^ 

B 4n 

+ ^(a)?-|f 2 exp( + iCi-iC 2 )_).] . 

The integration on and is taken over -ir/2 to tt/2 and, since 

% 'V ^ 

C-j and ? 2 are uniformly distributed, the last two averages are zero, 
so that: 

<V 2 > = C 2 [<g(0)? 2 >+^(«)f 2 >]. 
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Rewriting <v > in terms of the integrals and pdf ! s: 


<» 2 > 


? 2 

1 80 .a r .2 

2 L 
Es 4tt 


TT 

( , ? 

g(u)n“h^(n)<5 (a)dadn+ 


g(a)f,h (a)6( n )dadn 


Integration of the delta functions gives: 


<v 2 > = I|oV c 2 
B 4ir 


n 2 h n (n)dn+ \ g(a)f 2 h a ( a )dc 


Substituting the expressions for h (n), g(a), and h a (a). 


<v > 


180 

“IT 

B 4tt 



exp[-(n-f-| ) 2 /2a 2 ]dn 



TT 

^ exp 
o 


a 2 4&n2 

a 

2 

a 

t 

2 — 

2 eXp 
°2 

■ rr 

L 2o d 

da 


The first integral is the second moment of Gaussian variable ?-j , or 
<? 2> = f^(l+a 2 ), and the upper integration limit on a may be taken 
to so that: 


o2 
<v > 


180 2 A 2 r 2 
L 

B 4tt 


oa 

i (i+a i )+ il V exp 

i °2 


(8tn2a 2 +B 2 )a 2_1 
2a *BT 


da 



8tn2a 2 +B 2 

_ 

20^ 


Let 



For a Rayleigh distribution. 


1 f 

A 


r 2 

2 j a expL-a /2A“Jda 


A 2 = \ a exp[-a 2 /2A 2 ]da . 


Rewrite the integral in the expression for <v 


\ 2 . 


<v > 


2 2 

180 r 2 

"7 L 

B £ 4tt 


f 2 


f 2 (l+a 2)+ ( a exp[- 

°2 o 


The integral may then be replaced by A : 


<v^> = c 2 

B 4tt 


Tf(1+cr2) +T | h_ 
° 2 _ 


Finally, substituting in the definition of A , 


<^ 2 > = 1 80 2 A 2 c 2| 

B 2 4tt 


4 1 + ° 1> + 1 (— 4-2 

1 1 ^ \8An2aJ+B‘ 1 


in terms of A 


a 2 /2A 2 ]da 


one obtains: 




APPENDIX D 
MODEL SUMMARY 

I. Intensity Functions 

<? 2 > + <f^> = 1 

= f^O+^) 

f~ = 7? = 1 - exp[-L/L ]; L M80 km 
Id r o o 

L = path length (km) . 

II. Variance Parameters 

Amplitude Variance 

a 2 = 2.6 x 10" 7 f 7/12 (GHz) L 11/6 (km) . 

Angle of Arrival Variance 

a ^ = 5.67 x 10" 6 “^ 6 (km) d~ 1/3 (m) 

L = path length 

d = circular aperture diameter 
f = frequency 
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III. Gain Reduction 


< w> 

10 I og q — 

<v> 


0^=0 (no angle fluctuation) 


<v> = ensemble average of receiver voltage 


f? + 7 


, J 1 

1 T 2 


R = 10 log iq 


B" 

\4^n2a^+B^ 


—2 — 2 

T 1 T 2 


B = half power beamwidth in degrees assuming a Gaussian 
antenna pattern function. 


IV. Received Signal Variance 

2 ? 
c 2 in <v > - <v> 

S —10 log nn o 

10 < v > 2 


2 

5 0 i + 


—2 2 
f 2 B 


S 2 = 10 log 




f 2 

m " t 2 


B 


^4^n2o2+B 2 


10 


f 2 + f 2 
r l '2 


B 

4&n2a 2 +B 2 



APPENDIX E 

FADE DISTRIBUTION VARIANCE 


If 


S 2 = 10 log-jQ 


N 

I (V- - <v>) 

i=l 1 

N <v> 2 


is the definition of received signal variance, and we are interested in 

2 

the variance of the log amplitude, denoted o^, a transformation must be 
made. Let ^ be defined as log amplitude by 

= 20 1 °g -] 0 ^7 

If the fluctuations about <v> are small compared to the magnitude of 
<v>, let 


<v> 


v.-<v> 

= 1 + — = 1 +A ■ 

<v> i 


A. « 1. 


Then , 

A. = 20 log 1Q (1+A.) 

may be approximated by 

% (20 i °g i 0 e ) •A i . 

and 

Gaussian , 
log- 


If a function f(v) is normally distributed, with mean <v> 

O 

variance a the function kf(v), where k is a constant, is also 
v 2 2 

having mean k<v> and variance k a . From the definition of the 
3 v 

amp! itude. 
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v i =<v> 


*« • = k*. 

20 l og Q e i i 


However, the term on the left is a function of v. for which we have 

2* 1 2 

previously defined variance in dB as S . Hence, if is the variance 
of * i# 


S = 10 1 og-j 0 


20 log 1Q e I v t 


S = 20 log 


10 


( <*1 } 

20 lo 9]0 e , 


Solving for a , the standard deviation of the log-amplitude. 


a z = 20 log-| Q e • 10 


S 2 /20 


*Note: 


Var 


v .-<v> 

<v> 


<v> 


2 Var ^v i --<v> 


Var 


v . -<v> 

j 

<v> 


<v> 


2 " Var < v i 



for a symmetric density function. 


v.-<v> 


10 log 1Q Var \-~r S‘ 


V* 



